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We investigate the time evolution of optical conductivity in the half-filled one-dimensional ex¬ 
tended Hubbard model driven by a transient laser pulse, by using the time-dependent Lanczos 
method. Photoinduced in-gap excitations exhibit a qualitatively different structure in the spin- 
density wave (SDW) in comparison to the charge-density-wave (CDW) phase. In the SDW, the 
origin of a low-energy in-gap excitation is attributed to the even-odd parity of the photoexcited 
states, while in the CDW an in-gap state is due to confined photogenerated carriers. The signature 
of the in-gap excitations can be identified as a characteristic oscillation in the time evolution of 
physical quantities. 

PACS numbers: 71.10.Fd,03.67.Mn,75.10.Kt,75.25.Dk 


I. INTRODUCTION 


Time-resolved spectroscopy can provide new insights 
into the dynamical properties of condensed matter. 
Among many time-resolved spectroscopies, pump-probe 
optical measurements have been performed extensively to 
investigate strongly correlated electron systems. Pump¬ 
ing pulse applied to one-dimensional (ID) Mott insula¬ 
tors with spin-density-wave (SDW) correlation creates 
not only a photoinduced metallic state^i^ but also a pho¬ 
toinduced charge-density-wave (CDW) stated. In two di¬ 
mensions, the photoexcitation of underdoped cuprates 
induces a signature of enhanced superconductivity^^— . 
Even in a conventional superconductor, an observation 
of the amplitude mode^ is under intensive investigation 
within the field of nonequilibrium and nonlinear phenom¬ 
ena induced by a transient electric field. 

The time evolution of optical conductivity in the ultra¬ 
fast process has been theoretically examined^i^^— , since 
this quantity contains information on the photoinduced 
charge dynamics and the initial relaxation through inter¬ 
nal degrees of freedom of a given closed system. In Mott 
insulators, photoexcitation is expected to induce a new 
excitation channel inside the Mott gap. The nature of 
those states has not been fully understood from a theoret¬ 
ical viewpoint even for ID Mott insulators where, due to 
spin charge separation, the spin degrees of freedom affect 
the charge dynamics much less than in two-dimensional 
(2D) Mott insulatorsi^iii. 

In this paper we address the physical properties of pho¬ 
toinduced states, their positions within the Mott gap, 
and their influence on the time evolution of the nonequi¬ 
librium optical properties. For this purpose, we calculate 
the time-dependent optical conductivity of a ID extended 
Hubbard model driven by a transient laser pulse. In this 
model, one can change the ground state from the SDW 
phase to the CDW phase by changing the interaction 


parameters and tuning the laser pulse so as to enhance 
charge correlation even within the SDW phasei^. We 
find that photoinduced in-gap excitations appear in a 
different manner, depending on the nature of the initial 
ground state, i.e., SDW or CDW state. In the SDW, 
a dominant in-gap excitation emerges at the low-energy 
region reflecting the energy difference between the even- 
odd parity of photoexcited states. In the CDW, in-gap 
states appear in the middle of the Mott gap and their ori¬ 
gin can be attributed to photogenerated carriers. We also 
find that the signature of the in-gap excitations emerges 
in the time evolution of physical quantities such as the ki¬ 
netic energy and peak intensity of the in-gap excitations 
as a characteristic oscillation. 

The rest of the paper is organized as follows. In Sec. 
II, the model is introduced and the numerical method 
motivated by the experimental pump-probe setup is de¬ 
scribed. In Sec. Ill the characteristic optical responses 
of the system driven by a laser pulse are presented. The 
photoinduced in-gap excitations, which exhibit qualita¬ 
tively different structures in the SDW and CDW phases, 
are analyzed in detail. Their influence on the time evolu¬ 
tion of various physical observables is discussed. Finally, 
a summary is given in Sec. IV. 


II. MODEL AND METHOD 


We investigate the extended one-dimensional (ID) 
Hubbard model at half filling. 


H = -th'^ (cLci+i.o- + H.c.) 
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where cj ^ (ci,a) is the creation (annihilation) operator 
of an electron with spin cr at site i, Ui = rii^^ + 
is the hopping constant, and U and V are on-site and 
nearest-neighbor repulsions, respectively. 

The application of a spatially uniform and time- 
dependent electric field expressed by a vector poten¬ 
tial A{t) is incorporated into the hopping terms via the 
Peierls substitution: 

c\ ,^Ci+i^a + H.c. -b H.C., (2) 

leading to —?> H{t). We assume that A{t)^ representing 
the pumping pulse, is expressed in terms of the temporal 
gauge 
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^pump(t) = A^e cos [wo (t - to)], (3) 

where the temporal distribution of Apump(t) is Gaussian 
centered around tpi while the width is controlled by the 
parameter td, and wo is the central frequency. The aver¬ 
age incoming number of photons per lattice site during 
the pump is estimated to be oc ApWctd- 

In our calculation, the system’s temporal evolution 
is evaluated by the time-dependent Lanczos method, 
where, within the reach of the exact diagonalization, the 
time-dependent wave function tplt) can be obtained ex- 
actlyi^. Based on the information of '0(0, we calculate 
the time-dependent current density j{t) = (0(O|j|0(O), 
where j is the current operator. 

Motivated by the experimental pump-probe setup, we 
use the following strategy to obtain the time-dependent 
optical conductivity after pumping. For a given pump, 
we trace the temporal evolution of the system twice: 
first, without the probe and the second time with the 
applied probe. The current induced by the probing 
pulse Jpr(0 is obtained from the current difference be¬ 
tween these two consecutive simulations. The subsequent 
Fourier transformations on jpr(t) and Apr(f) produce 
the time-dependent optical conductivity^ as a {uj, At) ss 
jpr{z)/izAp^{z)^ where At is the time difference between 
the probing pulse and the pumping pulse and z = u! + iS 
with a small positive number S. Introduction of S is nec¬ 
essary for obtaining the spectral weight that is due to 
the Drude peak in the limit of w —>■ 0. We have checked 
that cr (u>, At) gives a good agreement with that obtained 
by a nonequilibrium linear-response theory!^ as long as 
the strength and the width of the probing pulse are small 
enough. We note that the result is also insensitive to the 
details of the probing pulse, regardless of the shape of 
the vector potential and the central frequency. 

In our numerical calculations, we consider a half-filled 
periodic lattice with L sites. We set = 1 as a unit of 
energy and h = 1. Throughout this work, the value of U 
is chosen to be t/ = 10, which is in the strong coupling 
regime. The central frequency of the pumping pulse wq is 
always set to match the absorption-peak energy of a given 
system. The small positive number S for the calculation 
of <t{uj, At) is taken to be <5 = 1/L. 


FIG. 1. (Color online) Re ( 7 ( 0 ;, At) for the SDW phase of the 
extended Hubbard model with U = 10. (a) At dependence for 
L = 14 and F = 3. Ao — 0.3 and td = 0.5. (b) To dependence 
for L = 14 and V — 4.5. At = 6 and td = 0.5. (c) L — 14, 
V = 4.5, Ao = 0.07, and At = 20. td = 5 is ten times larger 
than the case of (b), resulting in an enhancement of charge 
correlation after pumpingi^. (d) Comparison between L = 10 
and L = 14 with V = 4.5, To = 0.3, At = 6, and td ~ 0.5. 


III. RESULTS AND DISCUSSION 
A. In-gap excitations in SDW 

Let us start with the case where the initial ground 
state is the SDW state. Figure [TJa) shows the real part 
of the time-dependent optical conductivity, Recr(w, At), 
for a L = 14 ring with F = 3. Before pumping, there is a 
single absorption peak at w = 6.1, which is due to the ex- 
citon. With pumping, the exciton absorption peak loses 
its weight, while a low-energy peak at w = 0.2 emerges in 
the Mott gap. The low-energy peak persists throughout 
the time evolution. Figure [IJb) shows Tq dependence at 
At = 6 for F = 4.5. As expected, with increasing Tq, 
the exciton peak at a; = 4 loses its weight continuously 
and the low-energy in-gap peak increases. 

We now discuss the origin of the low-energy peak seen 
in Figs.IIJa) andITJb). Since the peak appears only after 
excitation with the pump pulse, its origin should be re¬ 
lated to photoexcited states. The photoexcitation excites 
carriers into the state that gives rise to the absorption 
peak in the equilibrium, which is an optically allowed 
odd-parity state. The probe pulse couples in part to the 
photoinduced carriers in the odd-parity state, resulting in 
an excitation from the optically allowed state to an op¬ 
tically forbidden even-parity state. In Mott insulators, 
it is known that the optically forbidden excitonic state 
is located slightly above the optically allowed excitonic 
state2£i2i. Therefore, nonequilibrium optical conductiv¬ 
ity should detect such an excitation. This is the origin 
of the low-energy peak. Its energy thus corresponds to 
the energy separation between the optically allowed state 
(exciton peak) and the forbidden state. In order to con- 
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FIG. 2. (Color online) (a) Contour plot of Re(T(aj,At) 
(1 < Ci) < 10) for the SDW phase of the L = 10 extended 
Hubbard model with [7 = 10 and V = 3. The amplitude is 
normalized by the maximum value of the equilibrium absorp¬ 
tion peak. The conditions of the pumping pulse are Aq = 0.3 
and td = 0.5. (b) The amplitudes of Fourier components of 
the maximum intensity (red dashed line) near cu = 6 and the 
corresponding position (black solid line) in (a). The ampli¬ 
tude is normalized by the maximum value. The blue thin line 
represents Re ( 7 ( 0 ;, At = 6) for the same system. 


firm this, we perform the same calculation on a smaller 
system size. The optically forbidden state is known to 
be lower in energy than the allowed one if the small sys¬ 
tem has periodic boundary conditions. A L = 10 lattice 
with V = 4.5 is such a case. If the optically forbidden 
state is lower in energy than the optically allowed one, a 
negative-energy excitation appears after pumping, lead¬ 
ing to negative weight at the positive-energy side since 
uja(uj) is odd in cu. This actually occurs in the L = 10 
lattice as shown in Fig. [Ifd), where Rea(uj,At = 6) is 
compared with that of the L = 14 lattice. The data for 
L = 10 exhibit a peak at w = 0.4 with a negative weight 
in contrast to the L = 14 case with a positive weight. 
This confirms our hypothesis that low-energy in-gap ex¬ 
citation comes from the excitation from the optically al¬ 
lowed to forbidden state. 

It is vital to test whether the photoinduced low- 
energy in-gap excitation demonstrated in Figs. lUa) and 
mb) can be detected by pump-probe optical experi¬ 
ments. A recent time-dependent optical reflectivity 
measurement on 6js(ethylendithyo)-tetrathiafulvalene- 
difluorotetracyanoquinodimethane ET-F 2 TCNQ has in¬ 


dicated the presence of photoinduced low-energy states 
located around ~100 meV by analyzing the intensity os¬ 
cillation near an exciton peab^. This energy is well below 
the energy of the exciton peak, 675 meV. The situation 
looks similar to Figs. [iKa) and (b). In order to explore 
the possibility of detecting the in-gap excitation, we cal¬ 
culate Re ( 7 ( 0 ;, At) and make a Fourier analysis, as is the 
case of the experiment. Figure UK a) shows the contour 
plot of Re < 7 ( 0 ;, At) (I < a; < 10) for a L = 10 ring 
with U = 10 and R = 3. We note that in this sys¬ 
tem and parameter set the photoinduced in-gap excita¬ 
tion emerges with a positive-spectral weight at w = 0.2. 
We can find an oscillating behavior of a photoinduced 
state below the exciton peak at a; = 6.3 as a function 
of At. The oscillation has a period of At 16. The 
Fourier amplitudes of the maximum intensity near uj = 6 
and the corresponding position are shown in Fig. HKb), 
together with Recr(a;, At = 6). We observe a dominant 
peak at w ~ 27r/16 ^ 0.4 in both the position and the 
intensity. The peak position of the Fourier amplitude 
roughly matches the position of the photoinduced in-gap 
excitation at a; = 0.2, which indicates it might come from 
the Rabi oscillation of the two odd-even states. We also 
find that such a low-frequency oscillation close to the en¬ 
ergy of the in-gap excitation does not occur in the cases 
with the negative-weight peak, e.g., R = 2 at A = 10 
(not shown). Therefore, we consider that the oscillation 
appearing in the exciton peak reflects the energy of the 
in-gap excitation. 


We next consider the case of a large pump width td = 5 
for V = 4.5, where a strong enhancement of short-range 
charge correlation after pumping has been reported in 
Ref. [ 1 ^. The large width of the pump pulse td leads to a 
well-tuned energy transfer between the pump pulse and 
the system. After the pump pulse we thus expect an en¬ 
hanced population of excited states around ojq, which is 
tuned to the position of the absorption peak. Such a non- 
thermal population of states around ujq naturally leads to 
the appearance of the negative weight of the absorption 
peak, as shown in Fig. [TKc). In addition to the negative 
weight, we find a pronounced structure around uj = 2 in 
Fig. [IKc) as compared to the case of small td [Fig. [IKb)] 
where an additional structure around oj = 1 exists. With 
increasing Aq up to 0.3 in Fig.lIJc), the structure around 
oj = 2 largely diminishes as it is accompanied by the 
vanishing of the induced charge correlations (not shown 
here). Since the pumping pulses with td = 0.5 used in 
Fig. HKb) also induce an enhancement of charge correla¬ 
tion but with less extent than that in Fig. [IKc)i^, the 
pronounced structure around a; = 2 in Fig. [TKc) may be 
attributed to a signature of enhanced charge correlation. 
In fact, the persistence of such an in-gap state is similar 
to the case of the CDW phase, as will be discussed below. 
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FIG. 3. (Color online) Recr(aj, At) for the CDW phase of the 
extended Hubbard model with U = 10. (a) At dependence 
for L = 14 and V — 6. Aq = 0.3 and td = 0.5. (b) Different 
parameter sets for L = 14 and H = 7. Aq = 0.3(1), td = 
1(0.1), and At = 20(2) for the solid (dashed) line. The in¬ 
gap peak with positive weight at cj = 5 is insensitive to the 
parameter sets, (c) Equilibrium optical conductivity Re ( 7 ( 0 ;) 
for the same parameter set as (a) but away from half filling 
(seven up-spin electrons and six down-spin electrons), (d) 
The same as (c) but the same parameter set as (b). 


B. In-gap states in CDW 

The CDW state is stable for V > 5. Figure[3Ka) shows 
the At dependence of Re ( 7 ( 0 ;, At) for R = 6. The weight 
of the absorption peak at w = 6.2 substantially decreases 
with pumping. A low-energy structure with both posi¬ 
tive and negative weights around ui = 0.3 emerges and 
exhibits a strong At dependence, which has not been seen 
in Fig. [H Because of a strong time dependence, we con¬ 
sider that the low-energy excitations are not a signature 
of photoinduced states. On the other hand, a strong in¬ 
gap peak appears at tu = 2.7 with weak At dependence. 
We regard this as a relevant photoinduced state. The 
Re ( 7 ( 0 ;, At) for V = 7 also shows a similar in-gap exci¬ 
tation that is independent of the pumping condition and 
At, as shown in Fig. EKb). Therefore, the in-gap state is 
expected as a characteristic of the CDW phase. 

We now discuss the origin of the in-gap state in the 
CDW phase. Again, we should notice that the pump¬ 
ing creates photo-carriers in the system. However, in the 
CDW phase, a metallic or low-energy excitation is un¬ 
stable as seen in Fig. |31 This is expected to be related 
to the charge distribution of the CDW ground state with 
the characteristic arrangement of alternating doubly oc¬ 
cupied and empty sites. The photo-carriers are created 
by transforming pairs of doubly occupied and empty sites 
to singly occupied ones, but the movement of electrons 
is still blocked by the presence of neighboring doubly 
occupied sites. This may lead to a finite-energy exci¬ 
tation inside the original CDW gap. This is confirmed 
by calculating the optical conductivity for the extended 
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FIG. 4. (Golor online) (a) At dependence of intensity of an 
in-gap state at tUin-gap = 4.45 in Re < 7 ( 01 , At) for the CDW 
phase of the L = 10 extended Hubbard model with U = 
10 and V = 7. The conditions of the pumping pulse are 
Ao = 0.3 and td = 0.5. (b) The amplitude of the Fourier 
component of the in-gap peak intensity (black solid line) as 
well as that of kinetic energy (red dashed line). The amplitude 
is normalized by the maximum value. The blue thin line 
represents Re <7(0;, At — 6). 


Hubbard model away from half filling with a single hole. 
Figures He) and [S^d) exhibit the equilibrium Ret7(a;) 
for U = 6 and U = 7 in the CDW phase, respectively, 
with seven up-spin electrons and six down-spin electrons 
in the L = 14 lattice. In-gap peaks appear at ui = 2.8 
and 4.6 for [7 = 6 and 7, respectively. Their energies 
are close to those of the photoinduced in-gap excitations 
in Figs. ISKa) and Hb), indicating that the presence of 
carriers is indeed the origin of the in-gap states. 

The intensity of the in-gap state at Win-gap = 4.45 for 
the L = 10 lattice with U = 10 and V = 7 oscillates with 
At as shown in Fig.HDa). After Fourier analysis, we find 
characteristic frequencies at w = 4 and 14 in Fig. Hb) 
(black solid line). The former is very close to the en¬ 
ergy of the in-gap state seen in Rea{uj,At = 6). We 
find that not only the intensity of the in-gap state but 
also the kinetic energy, which is an average value of the 
kinetic terms in H o, displays a pronounced weight in 
the Fourier spectrum around Win-gap, as seen in Fig. SDb) 
(red dashed line). Similar to the case of SDW, this char¬ 
acteristic oscillation related to the in-gap state is also 
recognized from the analysis of the remaining main ab¬ 
sorption peak (not shown here). This means that the 
energy of the in-gap state may be detected by measuring 
the oscillations of the physical quantities. The localiza¬ 
tion effect of the CDW background on the photoexcited 
charge carriers might be able to be identified as well in 
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higher-dimensional systems by investigating the super¬ 
imposed oscillations of time-resolved signals^^. The re¬ 
sponse at w = 14 is due to high-energy excitations that 
are due to symmetry restrictions not picked up by the 
optical conductivity measurement. These excitations are 
determined by U and V, leading to oscillations with a 
short period, as seen in Fig.|Tl[a). 

IV. CONCLUSIONS 

In summary, we have analyzed the physical properties 
of photoinduced excitations generated by a short laser 
pulse and their influence on the time evolution of various 
physical properties in the ID extended Hubbard model 
at half filling. We have calculated the time-dependent 
optical conductivity driven by a transient laser pulse and 
clarified the difference in the nonequilibrium optical re¬ 
sponse between the SDW and the CDW phases. The 
photoinduced in-gap excitations exhibit a qualitatively 
different structure in the SDW or CDW ground state. 
In the SDW, a dominant in-gap excitation emerges at 
the low-energy region, reflecting the energy difference of 
the even and odd parity of photoexcited states, where 
the even- (odd-) parity state corresponds to the optically 
forbidden (allowed) state. In the CDW, an in-gap state 
appears in the middle of the Mott gap and its origin 
can be attributed to photogenerated carriers. The sig¬ 


nature of the in-gap excitations in both cases emerges in 
the time evolution of physical quantities as a character¬ 
istic oscillation with a characteristic frequency roughly 
matching the energy of the in-gap excitations. Such a 
phenomenon has been reported in the time-dependent 
optical reflectivity measurement on ET-F 2 TCNQ, a pro¬ 
totype ID Mott insulator—. Observation of this charac¬ 
teristic frequency represents an experimental probe to 
identify photoinduced in-gap states. 
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